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ABSTRACT 



Off-center stellar tidal disruption flares have been suggested to be a powerful probe of recoiling 
supermassive black holes (SMBHs) out of galactic centers due to anisotropic gravitational wave 
radiations. However, off-center tidal flares can also be produced by SMBHs in merging galaxies. 
In this paper, wc computed the tidal flare rates by dual SMBHs in two merging galaxies before 
the SMBHs become self-gravitationally bounded. We employ an analytical model to calculate the 
tidal loss-cone feeding rates for both SMBHs, taking into account two-body relaxation of stars, 
tidal perturbations by the companion galaxy, and chaotic stellar orbits in triaxial gravitational 
potential. We show that for typical SMBHs with masses 10^ Mq, the loss-cone feeding rates 
are enhanced by mergers up to F ~ 10^^ yr^^, about two order of magnitude higher than those 
by single SMBHs in isolated galaxies and about four orders of magnitude higher than those by 
recoiling SMBHs. The enhancements are mainly due to tidal perturbations by the companion 
galaxy. We suggest that off-center tidal flares are overwhelmed by those from merging galaxies, 
making the identification of recoiling SMBHs challenging. Based on the calculated rates, we 
estimate the relative contributions of tidal fiarc events by single, binary, and dual SMBH systems 
during cosmic time. Our calculations show that the off-center tidal disruption flares by un-bound 
SMBHs in merging galaxies contribute a fraction comparable to that by single SMBHs in isolated 
galaxies. We conclude that off-center tidal disruptions are powerful tracers of the merging history 
of galaxies and SMBHs. 
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1. Introduction 

In the A-cold dark matter cosmology, both 
dark matter halos and galaxies form due to fre- 
quent mergers. In this paradigm, hierarchical 
galaxy mergers would incorporate multiple su- 
pcrmass ive black hol es (SMBHs) into a galaxy 
([Volont cri et al.ll2003[ ). When two SMBHs, ini- 
tially embedded in the two cores of the merging 
galaxies, sink to the common center of the system 
due to dynamical friction and become gravita- 
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tionally bound, a superm assive black hole binar y 
(SMBHB) would form (jBegelman et al.l Il980[) . 
During the interaction between the SMBHB and 
the stellar and gaseous environments, if the two 
SMBHs could successfully evolve to a separation 
of hundreds of Schwarzschild radius, then gravita- 
tional wave (GW) radiation could lead to the coa- 
lescence of the SMBHs within a Hubble time, and 
the asymmetry of GW radiation is predicted to im- 
part a recoiling velocity on the post- merger SMBH 
(|Hughesl[2009l : iGentrella et al.ll2010l ). Detection of 
the GW radiation from coalescing SMBHB would 
be a vital test of the theory of general relativity 
(GR), and is the major goal of the ongoing Pul- 
sar Timing Array (PTA) project and any future 
space-based GW mission. 

Despite many efforts to detect GW radia- 
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tion from coalescing SMBHBs, theoretical stud- 
ies found large uncertainties for the dynami- 
cal evolution of SMBHB in normal galaxies: in 
the absence of gas and efficient stellar relax- 
ation, the evolution of SMBHB would stall at 
sub-parsec (pc) s cale and not enter the GW 
radiation regime jMerritt fc Milosavlievil l2008t 
Colpi fc Dottill201l[ ). while recent jV-bodv simu- 
lations suggest that efficient repopulation of stars 
to t he galaxy core m a y be norm i n rea l merg- 



ers (jPreto et alJ I2OIII: Khan et all l201lh . Ob 



servationally, it is difficult to test the dynami- 
cal evolution of SMBHB in stellar systems, be- 
cause of the lack of electromagnetic (EM) radi- 
ation from the vicinity of the dormant SMBHs. 
Recently, "tidal flares", the EM outbursts pro- 
duced due to tidal disruption of stellar objects by 
SMBHs, have been identified as powerful probes 
of the mass and spin of the otherwise dormant 



SMBHs (Rees 19881 Komossa 


2002; Donley et al. 


2002; van Velzen et al. 2011; 


Gezari et al. 20091 


Bloom et al.l2011; Burrows et al.l2011; Cenko et al. 


2012). The flaring rate for 


a single SMBH in 



an isolate d galaxy is estimated to be 10 ^ to 

( Maeorrian fc Tremainell999t Sver fc Ulmer 



lO"'' yr"^ 



1999l:IWang fc Merrittl2004l:rBrockamp et al.l2011^ 

It is predicted that the formation and evolu- 
tion of bound SMBHB at galaxy center would 
significantly change the event rate and affect 
the light curves of tidal flares. Shortly after 
the formation of SMBHB, the three-body inter- 
action between the binary and a bound stellar 
cusp w ill enhance the flaring rate to as h i gh as 



(Ivanov et al.ll2005t IChen et al.ll2009l 



2011 



1 vr-i TV---. ' 

Wegg fc Nate Bodd 120111 ). After the stellar cusp 
is disrupted, mainly due to slingshot ejection, if 
stellar relaxation is inefficient in galaxy center, 
the flaring rate will become one order of mag- 
nitude l ower than that in the single black hole 



system (jChen et all l2008f ). When the SMBHB 



enters the GW-radiation regime, the interruption 
and recurrence of tidal-flare light curve by the 
perturbing seconda ry black hole o ccur on an ob- 
servable timescale ( Liu et a 1 120091 ). and the stars 
resonantly trapped by the inspiralling SMBHB 
may produ ce a tidal flare around the coalescence o f 
the binary (|Schnittmanll2010HSeto fc Mutoll201l[) . 
After the coalescence of the binary, the launch 
of a recoiling SMBH might be also accompanied 
by a brief burst of tidal flares. As the recoil- 



ing SMBH travels outside the galaxy core, tidal 
disruption of the stars gravitationally bound to 
the hole may produce a flare apparently displaced 



Bpiacea 

from the galaxy center (iKomossa fc MerrittI 20081; 
Merritt et al1l2009t IStone fc Loebll2012t iLi et al 



20121 ). Because of the many differences between 
the flaring rates in single and binary SMBH sys- 
tems, it is suggested that tidal flares can be uti- 
lized to constrain the fraction and dy namical evo- 



lution of SMBHBs in galaxy centers (iChen et al 



20081) . 

Off-nuclear tidal disruption flares and ultra 
compact star clusters with peculiar properties are 
suggested to be the key features of gravitational 
recoiling SMBHs in galaxies. However, off-center 
tidal disruption flares can also be produced by 
SMBHs in merging galaxies, and star clusters with 
the proposed peculiar characters may also form 
by tidal truncation of secondary galaxies in mi- 
nor mergers. In particular, tidal flaring rates 
would be enhanced at a stage when the SMBHs 
are still isolated in the cores of merging galax- 
ies because of the mutual perturbation between 
merging galax ies, ni u ch ear lier than the formation 
of SMBHBs. IRoosI (|l98l[ ) pioneered the discus- 
sions on the stellar tidal disruptions by assuming 
that merging galaxies harbor Sgr A*-like SMBHs 
and by taking into account perturbations by com- 
panion galaxies. However, it is unclear by how 
much the tidal disruption rates can be boosted in 
physical galaxy models combining the correlations 
of central SMBHs and galactic bulges, how many 
tidal flares in the universe are contributed by this 
merger phase, and how they would affect the con- 
straint on the merger history of SMBHs. As a 
first step toward addressing these issues, in this 
paper we calculate the stellar-disruption rates dur- 
ing galaxy mergers and investigate the prospect 
of using tidal flares to probe multiple SMBHs in 
merging galaxies. 

The outline of the paper is as follows. In § [21 
we introduce the basic loss-cone theory and the 
stellar-disruption process in single SMBH systems. 
In § [HI we describe the stellar relaxation pro- 
cess in merging systems and generalize the loss- 
cone theory to calculate the corresponding stellar- 
disruption rate. We also discuss our results for 
different merger parameters. Based on the cal- 
culated rates, we investigate the contribution of 
tidal fares by merging galaxies in § [4| and discuss 
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our results and their implications in § [5l 

2. Loss-cone feeding in single SMBH sys- 
tem 

We first calculate stellar-disruption rates for 
isolated galaxies with single SMBHs, to prepare 
the basics for more complicated calculations for 
merging galaxies. A star with mass m* and ra- 
dius r* would be tidally disrupted when it passes 
by an SMBH as close as the tidal radius 



n 



Ml 



1/3 



4.9 X 10""^ M- 



-1/3 



pc 



(1) 



(2) 



(|Hillsl I1975I: lReeslll988l) . where M, is the black 
hole mass, Mr = A/./lO"^ Mq, and Rq and Mq 
are, respectively, the solar radius and mass. In 
the following, we assume = Rq and m* = Mq 
unless mentioned otherwise. For these solar-type 
stars, when M, <C 4 x lO'^ Mq, tidal disrup- 
tion happens outside the marginally bound or- 
bit of the black hole, and collisions between the 
bound stellar debris, as well as the subsequent 
accretion on to the black hole, could produce a. n 
EM flare known as the "tidal flare" (|Reeslll988l) . 
The criterion for stellar disruption is then J < 
Jtd — (2G'M,rt)^/^, where G is Newtonian grav- 
itational constant, J is the specific angular mo- 
mentum, and Jtd is the specific angular momen- 
tum corresponding to a pericenter distance of r^. 
Here, the latter approximation accounts for the 
fact that most stars are disrupted along parabolic 
orbits, i.e., their specific binding energy £ <C 
GM,/rt. When M, > 4 x lO"^ Mq, the marginally 
bound orbit of Schwarzschild black hole becomes 
greater than the tidal radius, then the criterion 
for stellar depletion becomes J < Jmb, where 
J,nb denotes the specific angular momentum for 
marginally bound geodesic. In general, Jmb is a 
function of black-hole spin and inclination relative 
to the equatorial plane, but for simplicity we adopt 
the orientatio n-averaged value Jmb = 4:GM,/c 
(|Kesdenll2012l ) in the following calculation, where 
c is the speed of light. For even greater black- hole 
mass M, > 10^ Mq, tidal disruption occurs inside 
the event horizon of the central SMBH even when 



the black hole is maximally spinning, so no tidal 
flare could be produced by disrupting solar-typ e 
stars ( Ivanov fc ChernvakoTOll2006t lKesden|[20ia) . 



As a result of tidal disruption and direct cap- 
ture, a small fraction of stars are lost from the sys- 
tem during their pericenter passages. In a spher- 
ical system, the disruption rate of stars from dis- 
tance r to r + dr from the central SMBH is ap- 
proximately 



dT 



AT:r^drp{r) 0'^{r) 
TO* td{r) 



(3) 



(|Frank fc Reeslll976t [Sver fc Ulmeilll999l ). where 
p{r) is the stellar mass density at r, td{r) is the dy- 
namical timescale, and 0^(r) estimates the fraction 
of stars subjected to lose from the system. The loss 
fraction 0^ is dimensionless and can be interpreted 
geometrically as a solid angle, because at r the 
lost stars have velocity vectors pointing toward the 
SMBH within an angle of 9\c{r) = Jic/Jcir) and 
in an isotropic system their fraction is 9^ = 9f^. 
Here Jc denotes the angular momentum for cir- 
cular orbit and is of order r(7{r) given the stel- 
lar velocity dispersion a{r). The cone-like re- 
gion with half-opening angle 9ic toward the SMBH 
is therefore called "loss cone". The isotropy of 
stellar distribution breaks down at the edge of 
loss cone when the orbital-averaged rms veloc- 
ity deflection angle 9d(r) is much smaller than 
6>i. (iLightman &: Shapirol Il977t ICohn k Kulsrud 
Il978l ). Taking this effect into account, careful 
analysis of the loss-cone structure suggests that 



9' =mm{9l,9'jln9-^ 



(4) 



(|Young|ll977h . Therefore, when 9d ^ ^ic ( "pinhole 
regime"). Equation ^ recovers 9"^ = 9f^, because 
the stars act as if the loss cone does not exist and 
the system remains isotropic. On the other hand 
when 9d <^ 9\c ("diffusive regime"), the loss cone 
becomes empty within one dynamical timescale, 
so afterwards only a fraction ^^/jln^id of stars 
residing at the boundary layer 9\c ^ 9\c + 9d of the 
loss cone will be depleted during one td- The total 
stellar disruption rate F is an integration of Equa- 
tion ([3]) over both pinhole and diffusive regimes. 

To calculate p{r), td{r), and 9'^{r), a physical 
model describing the stellar distribution in the 
host galaxy needs to be specified. We consider 
only the bulge component of a galaxy because it is 
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the major source for stellar disruption. We model 
a bulge with a spherical model with double power 
laws, i.e., 



p(r) = 







(r < n) 

{n <r < r,nax) ,(5) 
(r > r,nax) 



where is the break radius, ph is the stellar 
mass density at r^, 7 and /3 arc, respectively, the 
inner and outer power-law indices, and r^ax is 
the cut off radius to prevent divergence of the 
total stellar mass. The five model parameters, 
('^fe, ''max, Pb, 7, are determined by the follow- 
ing five physically motivated conditions 

1. We define rf, as the influence radius of 
such that the enclosed stellar mass 

is 2 A/,. 

2. and 3. The values of 7 and /3 are adopted from 
empirical galaxy models (jFaber et al.lll997 : 



Lauer et al.l 120051 ) and will be specified ex- 



plicitly in the following calculations. In our 
fiducial model, 7 = 1.75 and /3 = 2, so that 
the galaxy has an inner B ahcall-Wolf and 
outer isothermal profile (jBahcall fc WoH 
19761 ). By varying 7 and /3 (7,/? < 3), our 



simplified galaxy model could reconcile with 
a variety of real galaxies. 

The total stellar mass enclosed in the ra- 
dius fmax is AM,, where A = 400 so that 
the SMBH-to-galaxy mass ratio satisfies the 
empi rical correlation in the local university 
(e.g.. lMarconi fc Huntll2003[ ). 



At the effective radius re, where the two- 
dimensional (2-D) surface-density isophote 
encloses half of the total galaxy mass, the 
stellar velocity dispersion satisfies the em- 
pirical correlation M, p i 10® ((Te/200 kms-^)-* M(i 
(jTremaine et al.ll2002h . Note that the stel- 
lar mass enclosed by the 3-D sphere of ra- 
dius re is Af^re) ~ (0.36, 0.32)^A/. when 
/3 = (2, 1.5), smaller than half of the galaxy 
mass. 



^It is suggested that when (/3 — 7) > 1, a pract ical definition 
for ri, is that the ma ss deficit inside ri, is 2M, l|Merrittl2006l : 

IStone fc Loebir2012l) . but the resulting pi, difli'ers from our 
fiducial value by only a factor of(/3 — 7)/(3 — 



According to Jeans's equation in the isotropic 
limit 



dipa^) Gp[M,{r)+M,] 



dr 



= 0, 



(6) 



the velocity dispersion a cx r^^/^ when r ^ ri, and 
a oc r^^^/"^ when r 3> r^; therefore, we calculate 
cr with 



cr(r) 



fT6(r/rfc)-i/2 < 
(Tb{T/nf'^'^ (r > n) 



(7) 



where Ub is the velocity dispersion at r^. By ap- 
plying Equations ([6]) and ([7]) at re, wc first derive 



^e + 1 GM, 

2/3-2 al ' 



(8) 



where A^ = M^{re)/M, and CTe is computed 
with Af, — (Te relation. Then the model pa- 
rameters (rh, rniax, Pb) are calibrated according to 
their definitions, and the results are rb ~ ''e[(6 — 

27)/(3Ae-/3^e)]^/(^-«, W (AMe)l/(3-/')re, 

and pb — (3— 7)A/,/(27rr^). For example, our fidu- 
cial galaxy model with Af, = 10^ A/q, 7 = 1.75, 
and /3 = 2 corresponds to rb ~ 4.5 pc, 7'e ~ 260 
pc, and 

' max 

~ 820 pc. 

Having specified the galaxy model, we now cal- 
culate the deflection angle 9d which determines 9^ 
in Equation ^ . Two-body scattering is an inher- 
ent relaxation mechanism in stellar system and it 
gives a lower limit of 62 = J2I Jc to 6*^, where Ji is 
the cumulative change of J due to two-body scat- 
tering during one dynamical timcscale. Because 
successive two-body scatterings are uncorrelated 
(incoherent), we have J2 = {td/trY^'^ Jc^ where 



t,(r) 



^/2a3(r) 



(9) 



nG^m^p^r) In A 
2V2B2 A/. 



(3 — 7) In A m„ 

-)7^)"^ (10) 

O-fc / \PbJ CFb 

is the two-body relaxation timescale. In A is the 
Coulomb logarithm (we assumed a fiducial value 
of 5), and 



B = 



3-7 



GM,/al (3 -/?)(/? -1) 



(11) 



is a correction factor of order unity. When two- 
body scattering dominates the relaxation process. 
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J2 is an increasing function of r, with the transi- 
tion between pinhole and diffusive regimes (J2 
Jic) being situated at r ^ rb. The differential loss 
rate dV/dr (eq. (O) scales as r^/^"^'' in the diffu- 
sive regime (r <^ rb) and as r~^~^/^ in the pinhole 
one (r 3> r^); therefore, the stellar disruption rate 
peaks at the transition regime at r ~ rf,. 

Take our fiducial model with M, ~ 10^ Mq, 
7 = 1.75, and /3 = 2 for example. The criti- 
cal radius where 62 = Of^ is rcri — 2.3rt, and 
the total disruption rate due to two-body relax- 
ation is r ~ 2.3 X 10~^ yr~^ , consistent with pre- 



vious calculations (e.g., iMaeorrian fc Tremaim 



1999"; 'Sv er fc UhTieilil999l : IWang fc MerrittI 



200 



Brockamp et al.l 20 111) . If M, increases, Tcn/ri, 



will also increase, given the fact that 9^^ is a dc- 



cx 



M: 



creasing function of r/r^, and that 
and 9^^ cx mI^^ at any r/r^. On the other 
hand, the integrated stellar-disruption rate will 
decrease, mainly because the diffusive regime of 
loss cone becomes larger. A more accurate cal- 
culation of r could be carried out by solving the 
diffusion equation in the 2-D £ — J space (e.g. , 
Lightman fc Shapirolll977l:ICohn fc Kulsrudlll978l: 



Magorrian fc Tremaind 119991 : IWang fc Merritt 
20041 ) ■ but it is considerably time-consuming and 
out of the scope of this paper. Nevertheless, the 
present scheme gives good approximation to the 
two-body disruption rate, and is sufficient to pro- 
vide references for the sake of investigating the 
effects of galaxy mergers on the stellar-disruption 
rate. 



3. Enhanced loss-cone feeding during galaxy 
merger 

Because the loss cone is already "full" in the 
pinhole regime, enhancing relaxation efficiency in 
this regime does not increase the fraction of loss- 
cone stars, therefore would not increase stellar- 
disruption rate. On the other hand, the loss cone 
in the diffusive regime is largely empty, so the 
disruption rate can be enhanced if stellar relax- 
ation in this regime becomes more efficient. En- 
hancement of stellar relaxation in the diffusive 
regime can be achieved by galaxy merger due to 
at least two processes. First, perturbation by 
the companion galaxy wou ld secularly change the 
stellar angular momenta (|Roodll98l[ ). Second, 
the triaxial gravitational potential built up during 



merger (jPreto et al.l201lHKhan et al.ll201l[ ) would 
drive stars to galaxy c enter in a chaotic manner 
(jPoon fc MerrittI I2OO II ). In this section we calcu- 
late the stellar-disruption rates due to the above 
two processes, and we show the rate for each of 
the two SMBHs in the merging system. 

3.1. Basic Theory 

A companion galaxy would tidally torque the 
stellar orbits in the central galaxy, secularly chang- 
ing the orbital elements. Given mass Mp of the 
perturber and its distance d from central galaxy, 
one can derive GMpv/d'^ for the tidal force exerted 
by Mp across a stellar orbit of radius r <C d in the 
central galaxy. The corresponding tidal torque on 
the stellar orbit is of magnitude Tp ^ GMpr^ / d^ . 
Because of the tidal torque, the angular momen- 
tum of star changes coherently, i.e., A J oc t, up 
to a timescale C where C is determined by the 
shorter one between the dynamical timescale of 
the perturber and the apsidal precession time scale 
of the stellar orbit ( Binnev fc Tremainell2008 ) . For 
t > ti^, the torque on stellar orbit adds up stochas- 
tically and in this case AJ^ oc t. Therefore, av- 
eraged over a timescale much longer than tuj, the 
tidal torque changes by an amount of Jp = 
Tptujtd{r) during each stellar dynamical timescale. 
As a result, the deflection angle 0^ in Equation (|4]) 
increases by an amount of 6'^ = (Jp/Jc)^- We 
note that the calculation of Jp is analogous to the 
calculation of angular-momentum change due to 
resonant relaxation where the resonance torque 
i s induced by the grainy gravitational potential 
(iRauch fc Tremaind Il996l : iHopman fc Alexander 
20061) . 



Galaxy merger also increa ses the triaxiality 
of the gravitatio n al potential (IPreto et al. 201 ll 
Khan et ahlboill ). IPoon fc MerrittI (|200ll ) showed 
that when the triaxiality is large, a consistent frac- 
tion of stars are fed to the loss cone in a chaotic 
manner and the loss cone remains full. Suppose 
fc is the fraction of stars on chaotic orbits, the ex- 
tra contribution to stellar-disruption rate can be 
calculated by replacin g 0^ with 9l = fc9y^ In 9\c 
(jMerritt fc Poonll20o3 ). It has been shown that 



fc approaches unity when the triaxiality becomes 
greater than 0.25, but will rapidly decrease to 
inside the influence radius of the central SMBH 
wher e the gravitational po tential is largely spher- 
ical (|Poon fc Merrittll2004l) . 
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Because of tidal perturbation and triaxiality 
during galaxy merger, the effective deflection an- 
gle 9'^ increases to 



(12) 



and in the diffusive regime the loss-cone-limited 
deflection angle (0^ in Equation ([4])) also becomes 
larger. Consequently, an enhancement of stellar- 
disruption rate is anticipated. Now we have pre- 
pared Equations ([3]), dl]), and ([T^ to calculate the 
stellar-disruption rate in merging galaxies. How- 
ever, the equations are valid only in the adia- 
batic approximation, i.e., the gravitational poten- 
tial varies on a timescale much longer than the 
typical timescale for stellar orbital evolution. If 
the adiabatic condition is violated, the galaxy 
core will be subject to significant hea ting and 
expansion on the dynamical timescale (lOstriker 



1972[lMerritt fc Cru j|200lt iBovlan-Kolchin fc Ma 
20071 ). For the stars at r ~ which predom- 
inate the loss-rate enhancement, the maximum 
timescale of coherent angular-momentum change, 
tijj, is limited by the apsidal precession timescale, 
which is of order td{rb)- The timescale for chaotic 
orbital evolution is also of order td{rb). The adi- 
abatic limit therefore requires the orbital period 
of the merging galaxies to be longer than tdiri,). 
For this reason, the following calculations are re- 
stricted to d > 2rfc. 

3.2. Stellar-disruption Rates 

We now calculate the stellar disruption rate for 
both galaxies in a merger. The black-hole and 
bulge components are modeled with the parame- 
ters (M,,7,/3), as is described in § [21 The mass 
ratio of the galaxies, by construction, equals the 
mass ratio of the SMBHs, q = M,,s/M,^rn < 1, 
where the subscript m denotes the quantity for 
the bigger main galaxy and s for the smaller satel- 
lite galaxy. As we have shown that the contribu- 
tion to stellar tidal disruptions is dominated by 
the stars at the break radius of galaxy, we can ap- 
proximately construct a merger system of galaxies 
without loss of generality as follows. Given the 
distance d between the two galaxy centers, the to- 
tal stellar density at any location is approximated 
by summing the densities of the two unperturbed 
bulges. In this density field, each galaxy approx- 
imately preserves its initial structure out to a ra- 
dius min(rmax, ''tr), where rtr is the truncation ra- 



dius due to mutual tidal interaction, defined by 
the condition that the mean densities within rtj- 
are the same for the two truncated galaxies. Fig- 
ure [1] shows the density contours (upper panel) , 
as well as the density distribution along the line 
connecting the two black holes (lower panel) , for a 
merging system with il/, = 10^ Mq, q = 0.3, and 
d = SOrfc. 

Given the configuration of the merging system, 
we calculated 62 and F due to two-body relaxation 
for each of the two galaxies. To calculate Op and 
F due to tidal perturbation, the perturber mass 
Mp is derived by integrating the stellar and black- 
hole masses in the perturber galaxy enclosed by 
rtr. Note that the perturber is the satellite galaxy 
when calculating F for the main galaxy, but can 
also be the main galaxy when calculating F for 
the satellite. To calculate 6c and F due to chaotic 
loss-cone feeding, the triaxiality of galaxy needs 
to be determined. But our model is axisymmet- 
ric by construction, so the triaxiality cannot be 
derived self-consistently. We circumvent this in- 
consistency by assuming that at any radius where 
the density increment induced by the perturber 
excesses 5 = 20% of the initially unperturbed den- 
sity, a fraction of fc = 50% of stellar orbits are 
chaotic. Otherwise /c = if 5 < 20%. The ra- 
dial range where fc = 50% is insensitive to the 
choice of S because of the steep density profiles we 
adopted in the following calculations. 

Figure [2] shows the stellar disruption rates as a 
function of d for both main (upper panel) and 
satellite (lower panel) galaxies. The parame- 
ters are {MT,q,-f,l3) = (1,0.3,1.75,2) by default. 
When d ^ lOOrf, « 450 pc, the loss-cone filling 
in both galaxies is dominated by two-body re- 
laxation (dotted lines) and the disruption rate is 
identical to that for isolated single SMBH. As the 
distance shrinks to d ^ lOOr^, about 2re of the 
central galaxy, the disruption rates induced by 
companion galaxies start to exceed those due to 
two-body relaxation. This is because 9p(ri,) be- 
comes greater than 6*2 (?'{,). As d further decreases 
to d < lOr;, « 45 pc, 9c{rb) becomes greater than 
(^2{'fb), so the contribution to F due to triaxial 
potential starts to exceed that due to two-body 
relaxation. When the two galaxy cores are as 
close as the break radius of the main galaxy, F 
in both galaxies have been enhanced by two or- 
ders of magnitude. In the subsequent evolution 
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Fig. 1. — Upper: density contour in the mid-plane 
of a merging system with A/, = 10^ Mq and 
q ~ 0.3. The two galaxies, both have 7 = 1.75 
and /3 = 2, are separated by 50rf, where rf, refers 
to the break radius of the main (bigger) galaxy. 
The dashed circles mark the tidal truncation radii. 
Lower: density distribution along the line connect- 
ing the two black holes (solid curve). The dotted 
lines show to the initially unperturbed density dis- 
tributions. 



with d < 2rb for which our simple scheme cannot 
be applied, the three-body interactions between 
the two gravitationally bound SMBHs and the 
surrounding stars are expected to play an impor- 
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In a real merger, becau se galaxy orbitals are ec- 
centric ( Jiang et al.ll2008[) . the distance d will not 
decrease monotonically, but oscillate between the 
apocenter distance rapo and pericenter distance 
Tper, both distances decreasing with time due to 
dynamical friction. In this case, one can average 
the stellar-disruption rate over one orbital period 
according to f = C'''"- T{r)v-^dr / f""" ^-^dr, 
where Vr denotes the radial velocity of galaxy at 
distance r. Since in our model, where r/vr oc r^/"^ 
and r(r) cx r~'', both /3 and 77 are of order 
unity (see Section [33]), we find that F differs from 
T{d) by a factor of also order unity if we define 
d = (rper + rapo)/2. In this sense, the rates in Fig- 
ure [2] can be used as the orbital-averaged stellar- 
disruption rates for galaxy mergers with eccentric 
orbits. We also note that we may have underes- 
timated the contribution from triaxial potential, 
because in our model by construction /c vanishes 
inside about the influence radius of black hole, as 
the stellar-density variation 5 inside the sphere of 
radius ri, is small (e.g. lower panel of Figure [T}. 
In real galaxies, however, chaotic orbits may par- 
tially exist inside the in fluence radius of black hole 
(IPoon fc Merrittll200l[ ). 



Dependence of disruption rate on 
model parameters 



3.3. 



In § 13.21 we have shown that tidal perturba- 
tion by the companion galaxy dominates the en- 
hancement of r in a merger. The enhancement 
occurs when Op in the diffusive regime exceeds 62- 
As a result, the critical radius r'cri,p that sepa- 
rates the pinhole and diffusive regimes is now de- 
termined by 0p(rcri,p) = ^ic, and enhancement of 
stellar-disruption rate requires that rcri,p < ^cri. 
Now we investigate in what mergers the condition 
''cri,p < ''cri would be Satisfied. 
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According to Jp(rcri,p) = Jic and the relation 

3(3-^ip)//3p 



— oc „3(l-/3p)/(2/3p) f ^ 



nj 



/3(/3p-3)//3p 



where 



C 



n 



B 



Ae(3-/3) 

6-27 



(2-/3)/(3- 



(13) 



(14) 



we first derive the following scaling relation in the 



limit 
galaxy: 



crijP 



< riy and J\c = Jtd for the central 



^ oc Sl/7c-l/7^,^-l/42g3(l-/Jp)/(7/3p) ^ 



6(3-/3p)/(7/}p) 



7-6 



2/3(/3p-3)/(7/3p)+6/7 

(15) 



Since 1 < (3 < 3 for the majority of galaxies 
( Lauer et al 1 120051) ■ Equation pj]) suggests that in 
general enhancement of F would occur when the 
perturbing galaxy is larger or the galaxy distance 
is smaller. When tidal perturbation dominates 
the loss-cone filling, according to Equation ^ and 
M, (X erf, the rate F in the limit rcn.p ^ rt, scales 
as 

F oc C-^/'B-'/Hi:/'\r,,,Jny/'-''. (16) 

For our fiducial model with 7 = 7/4 and /? = 2, 
we can derive F oc (d/rt,)"^/^, which is consistent 
with the numerical results given by the dashed 
lines at d < SOr^ in Figure [H 

Figure [3] shows the dependence of F on black 
hole mass when q = 0.3. Equations ([T5l) and ((T6)) 
suggest that F oc M,^'^'^'^^^'^^ when q and d/ri, 
are fixed. The enhanced stellar-disruption rates 
in Figure |3] generally agree with this scaling when 
M. < 4 X 10^ Mq. When M. > 4 x 10^ Mq, the 
dependence of F on M, steepens because direct 
capture of stars by SMBH (GR effect) becomes 
important, such that the scaling of loss-cone size 



changes from Jj^ oc M, 



4/3 



to Jl OC M2. When 
^» ^ 10** Mq, the loss-cone stars will be directly 
captured by the central SMBH without produc- 
ing tidal flar es if the SMBH is n on-rotating or ro- 
tates slowly ( Ivanov fc Chernvak ova 2006: Kcsdc3 
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Fig. 2. — Stellar disruption rates as a function of 
galaxy separation for main (upper) and satellite 
(lower) galaxies. The dotted, dashed, and dash- 
dotted lines refer to rates induced by, respectively, 
two-body relaxation, tidal perturbation, and tri- 
axial gravitational potential. The thin dashed 
lines show the analytical solution F oc dr^/"^ in 
arbitrary units derived in Section 13.31 The model 
parameters are {MT,q,'-f, (3) = (1,0.3,1.75,2) and 
refers to the break radius of the main galaxy. 
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Fig. 3. — Total steUar disruption rates as a func- 
tion of galaxy separation for different black hole 
masses. The model parameters are the same as 
in Figure El except that Mr = (100,10,1,0.1) 
from top to bottom with decreasing line thickness. 
The dashed lines indicate that central SMBHs are 
more massive than 10* Mq and stars fall into the 
SMBHs without tidal disruption. 



2012I ). and the corresponding curves are shown in 
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dashed lines. Note that even when the SMBH in 
the main galaxy is more massive than lO'^ Mq, the 
merging system could still produce tidal flares, due 
to the existence of a smaller SMBH in the satel- 
lite galaxy. We found that when 1 < M7 < 10, 
the event rates of tidal flares can be as high as 
^ 10^^ yr^^ as d shrinks to about rf,. 

Figure m shows the dependence of F on the mass 
ratio q = M,^s/M,^rn < 1 of the two black holes, 
while M,_„i is fixed. For both main and satellite 
galaxies, the enhancement of F becomes more sig- 
nificant as q increases. It is worth noting that 
even q is as small as 0.01, the stellar disruption 
rate in the main galaxy can still be enhanced by 
two orders of magnitude when d shrinks to about 
Tfe. We also find that the enhanced stellar dis- 
ruption rate in the satellite is more sensitive to q 
than that in the main galaxy. This is because the 
baseline stellar-disruption rate, i.e., the rate for 
single black hole in isolated galaxy, changes with 
q for satellite galaxy, but does not vary for the 
main galaxy since in the calculation Af,_,„ is fixed 
(e.g. see Equation [E]). Quantitatively speak- 
ing, according to Equations (fTS]) and (fT6|). when 
varying q while keeping M,,m fixed, the enhanced 
stellar disruption rate for the main galaxy scales as 
g3{i-PsKi-2-/)/{iiM^^]j[\c the rate for the satellite 

scales as g(24+7e)/42-(l-27.)(3+/3/3.-3/J.)/(14/^)^ p^^ 

example, given (7, /3, 7^, /J^) = (1.75,2,1.75,2), 
one can derive F oc q^^/^^ for the main galaxy 
and F cx gS^/s^ for the satellite. 

Figure [5] shows the variation of stellar disrup- 
tion rate when the density profile of the main 
galaxy changes. For the main galaxy, when the 
inner power-law index 7 decreases from 1.75 to 1, 
the stellar disruption rates due to two-body re- 
laxation and tidal perturbation both drop by a 
factor of a few, because of the slight decrement 
of the stellar density at r ^ rf,. Meanwhile, the 
dependence of F on d/rj, at d < lOr;, changes 
from (d/rfc)"^/'' to (d/rb)"^/^, resulting in an even 
smaller rate at d ^ 7't, . When the outer power law 
index (3 decreases from 2 to 1.5, the stellar disrup- 
tion rates in the main galaxies drop approximately 
by a factor of 20. This is because the galaxy 
with shallower outer density profile is more spa- 
tially extended and has lower central density. For 
the satellite, when 7 or /3 of the main galaxy de- 
crease, the enhancement of stellar disruption rate 
occurs at smaller d/rb and becomes weaker for a 
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Fig. 4. — Stellar disruption rates as a function of 
galaxy separation for different q. The other pa- 
rameters are the same as in Figure [H 



10"' 

10""^ 

10-= 

10"* 
10"^ 

18- 



^ 10" 



10- 





(7,/?) = (1.75,2| ^ 






r Main galaxy ~ 

; — ^ 




; Sat"elrit-e 







10 



100 



1000 



Fig. 5. — Stellar disruption rates as a function of 
galaxy separation for different density profiles in 
the main galaxy. The other parameters are the 
same as in Figured 
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fixed d/rfc. This is because of the main galaxy 
becomes greater as 7 or /? decrease, so that for 
the satclHte the physical distance of the perturber 
increases if d/r^ is fixed. We notice that when 
P = 1.5, the stellar disruption rate in the main 
galaxy remains lower than that in the satellite as 
d decrease. This result implies that in mergers 
where the main galaxies have low surface bright- 
ness, the tidal flares are mostly contributed by the 
satellite galaxies. 

When the density profile of the satellite galaxy 
is varying, the resulting stellar disruption rates are 
shown in Figure El In general, the dependence of 
r on the density profile can be understood in the 
light of the analysis for Figure El except that now 
the role between the main and satellite galaxies 
switches. However, one difference is that when 
d shrinks to about r;,, the disruption rate in the 
main galaxy is not sensitive to the density profile 
of the satellite. This is because when d ^ ri, the 
stellar cusp surrounding the SMBH in the satellite 
is almost completely striped off by the tidal filed 
of the main galaxy, so for the main galaxy the 
perturbing mass is approximately M, g. 

Figures j^HHl showed that galaxy merger starts 
enhancing stellar-disruption rate when the two 
galactic nuclei are still widely apart, well before 
the two SMBHs become gravitationally bound. 
The boost factor for each SMBH incorporated 
is about 102(M./10^ MQ){d/rby{q/O.^Y, where 
M, is the mass of the subject black hole, rf, 
refers to the break radius of the more massive 
galaxy, and fi and v are indices depending on the 
density profiles of the two galaxies. Less mas- 
sive black holes have smaller boost factors be- 
cause prior to merger they already have higher 
stellar-disruption rates. The exact boost fac- 
tor depends on the stellar-disruption rate prior 
to galaxy merger, which deserves some discus- 
sion. When calculating F for isolated galaxies, 
we considered only two-body relaxation but not 
more efficient relaxation processes, such as res- 
onant relaxation, perturbation by massive ob- 
jects, or relaxation pr ocesses in triaxial gravi- 
tational potential (e.g. Ranch fc Tremaind 1996t 
Perets et all 120071 : iMerritt fc PoonI \2004} . Res- 
onant relaxation enhances stellar-disruption rate 
only mildly, less than a factor of a few in typi- 
cal galaxies ( Ranch &: Ingalls 19981 ). Massive per- 
turbs, such as molecular clouds and stellar-mass 



black holes, if highly concentrate inside the influ- 
ence radius of an SMBH, in principle could en- 
hance t he stellar-disruptio n rate by orders of mag- 
nitude (jPerets et al.ll2007l ). But such galactic nu- 
clei could only be transient, because large concen- 
tration of massive perturbs normally corresponds 
to short relaxation timescale. On the other hand, 
weak triaxiality seems intrinsic to galaxies, sug- 
gesting that chaotic loss-cone feeding may be im- 
portant prior to galaxy mergers. If we use Equa- 
tion (119) in iMerritt fc VasilievI (120111 1^ to esti- 
mate the stellar-disruption rate induced by tri- 
axial potential inside the black-hole influence ra- 
dius, meanwhile use formulae derived in Sections^ 
and 13.11 with fc ~ 0.1 to calculate the rate due 
to chaotic orbits outside the black-hole influence 
radius, then the total disruption rates for iso- 
lated fiducial galaxies become F ~ (5.8,8.9,47) x 
10-5 yr-i when M, = (10*^,10^10*) Mq. For 
comparison, the rates due to two-body relaxation 
only are (4.0, 2.3, 1.4) x 10"^ yr~^. The difference 
is the greatest in the case of M, = 10* Mq, be- 
cause the "gap" between rcri and ri, is the largest. 
These results suggest that only in the most mas- 
sive galaxies with M, > 10* Mq could intrin- 
sic triaxiality make the enhancement of stellar- 
disruption rate less significant. 

4. Contributions of tidal flares by merging 
galaxies 

In a synoptic sky survey, the probability of 
catching tidal flares in merging galaxies does not 
depend only on the stellar-disruption rate, but also 
on the duration of galaxy mergers. In other words, 
the fractions of tidal flares in merging and in nor- 
mal galaxies are proportional to the numbers of 
stellar-disruption events produced during, respec- 
tively, the merger and the quiescent phases. Since 
the duration of a galaxy merger is determined by 
the dynamical friction timescale, tdt, the fraction 
of tidal flares in merging galaxies is proportional 
to the typical number of tidal stellar disruptions, 
n = tdfP. Given the distance d between two merg- 
ing galaxies, we calculate the dynamical friction 



iMerritt fc VasilievI 1 I2OIII) used 7 = 1.5 to derive the stellar- 
disruption rates. To derive rates for different 7, we used the 
scaling relation between r^rit ^'iid 7 above Equation (116) 
in their paper. 
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timescale as 



Ms HMg/Ms 



(17) 



see eq. [8.13] in Binnev fc Tremain3l2008[) . where 



Mg{d) here refers to the stellar mass enclosed by 
the radius d in the main galaxy and Ms is the 
total mass of the truncated satellite. When the 
two galaxies are distant and T is not enhanced, 
the total number of disrupted stars is proportion 
to the dynamical friction timescale, which is 



tdf 



q-hd{d) I HMg/Ms) (18) 



We refer to this early evolutionary stage as phase 
I. During phase I, the dependence of t^fV on d'^/^ 
implies that the majority of tidal flares are con- 
tributed by wide galaxy pairs. When loss-cone 
feeding is enhanced due to tidal perturbation by 
the companion galaxy, the stellar-disruption rate 
r increases with decreasing d. We refer to this 
later evolutionary stage as phase II, and we note 
that main and satellite galaxies enter phase II at 
different times. During phase II, the dependence 
of idfT on d flattens compared to that in phase I, 
implying a enhanced detection rate of tidal flares 
in close galaxy pairs. 

Figure [7] gives the typical number of disrupted 
stars (<dfr) function of d in our fiducial model. 
In the calculation, we did not consider the de- 
crease of stellar density due to tidal disruption, 
because the total mass of disrupted stars is neg- 
ligible with respect to the stellar mass in the ini- 
tial condition. In general, when the merger is in 
phase I, ttifT scales as d, as predicted. During this 
phase, tdfr is not sensitive to the total mass of 
the system as long as q is fixed, because massive 
systems where V is larger have shorter tdf- Note 
that before the main galaxy (solid curve) enters 
phase II, the satellite galaxy (dashed curve) con- 
tributes comparable number of, if not more, tidal 
flares. This is because when two-body relaxation 
predominates the loss-conc filling, smaller galax- 
ies have smaller diffusive loss cones, therefore will 
have higher stellar-disruption rates, as is explained 
in the end of Section [2l When the galaxy mergers 
enter phase II, which is marked by the dots, the 
fdfl" curve flattens, indicating an enhanced contri- 
bution of tidal flares by closer galaxy pairs. Dur- 
ing this phase, the contribution of tidal flares from 
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Fig. 6. — Same as Figure [Sj but varying (7s, /3s) 
in the satellite galaxy. 
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Fig. 7. — Typical number of disrupted stars con- 
tributed by main (solid) and satellite (dashed) 
galaxies at different separations. Lines with de- 
creasing thickness refer to systems with decreas- 
ing M7. The other parameters are the same as in 
Figure [2] 
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Fig. 8. — Same as Figure [7] but varying one model parameter, which is indicated at the upper-left corner of 
each panel. 



main galaxy is typically greater than that from 
satellite. 

Figure [5] shows the dependence of t^fT on differ- 
ent parameters of galaxy merger, which arc sum- 
marized as follows. 

1. Comparing panels (a) and (b), one can see 
that during phase I, t^fT scales as a 
characteristic relation due to dynamical fric- 
tion. In phase II, as q decreases from 1 to 
0.01, tdfF for the main galaxy increases by 
a factor of 10 relative to that in the fiducial 
case in Figure 7, but that for the satellite 
does not significantly change. 

2. When 7 of the main galaxies decreases from 
1.75 to 1 as shown in panel (c), in phase 
I the main galaxy contributes slightly less 
tidal flares compared to that in the fiducial 
model, because F is smaller as 7 decreases. 
Meanwhile, the satellite contributes slightly 
more tidal flares because of longer fdf- As 
a result, the relative contribution of tidal 
flares by the satellite becomes greater dur- 
ing phase I. During phase II, tdfT for both 
main and satellite galaxies increases as d de- 
creases from about lOOr?, to lOr;,. As d be- 
comes smaller than lOr;,, the typical number 
of stellar disruption decreases more steeply 
with smaller d compared to that in the case 
with 7 = 1.75, because of less enhancement 
of stellar disruption rate as shown in Fig- 
ure [3 



3. Panel (d) shows that when f3 of the main 
galaxy decreases from 2 to 1.5, the number 
of tidal flares from the main galaxy drops 
during both phases I and II by about one or- 
der of magnitude relative to that in the fidu- 
cial case, because F decreases significantly. 
While that from the satellite increases in 
phase I due to longer tdf and significantly 
drops during phase II because of greater 
physical distance between the galaxies. As a 
result, tidal flares arc dominantly from the 
satellite galaxies during phase I, and almost 
equally contributed by the main and satellite 
galaxies during phase II. 

4. Panel (c) indicates that when 7s of the satel- 
lite galaxy decreases from 1.75 to 1, the tidal 
flares contributed by the satellite become 
slightly less than those in the flducial case in 
both phases I and II. During phase II, when 
^''b ^ d < lOrf,, the number of tidal flares 
contributed by the main galaxy slightly in- 
creases relative to that in the fiducial case, 
because the satellite is more severely tidally 
truncated so that tdf becomes longer. 

5. When /3s of the satellite galaxy decreases 
from 2 to 1.5 as shown in panel (f), the 
tidal flares contributed by the satellite dur- 
ing phase I are one order of magnitude less 
than those in the fiducial case. During phase 
II, the main galaxy contributes more tidal 
flares than in the flducial case, because now 
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the satellite is more susceptible to tidal strip- 
ping and tdf is much longer than that in the 
fiducial case. 

Figures [7] and [8] suggest that during phase I the 
total number of disrupted stars scales roughly as 



ni{q,d) - I0^q-\d/10^nf/^ 



(19) 



insensitive to the total mass of the system or the 
stellar density profiles of the merging galaxies. 
During phase II, when the curve of t^fV is nearly 
independent of distance d, the total number be- 
comes about 



nii(g,A/7) ^ 200 q"'M^, 



(20) 



where m ~ 3(1 - /3s)(l - 2j)/{U/3s) - 1 and 
n = (45-|-7)/42 are power-law indices derived from 
Equations ([T5|). and ([TS)) . It is worth noting 
that each merger investigated above involves only 
two galaxies. However, mergers of group galax- 
ies are also common and the tidal disruption rates 
are expected to be even more heavily enhanced 
because of stronger perturbations and larger tri- 
axiality. 

5. Discussions 

Formation of SMBHBs at galaxy centers is an- 
ticipated in the paradigm of hier archical galaxy 
formation ( Begelman et al.l llQSOl ). and coales- 
cence of the binaries is predicted to induce re- 
coiling velocities on the post-merger SMBHs 



Centrella et"all l2010l). In ou r previous works 



Chen et all l2008l l2009l I2OIII: iLiu et all l2009l) . 

investigated the possibility of using tidal- 



we 



disruption flares to identify gravitationally bound 
SMBHBs of sub-pc separations in galactic nuclei. 
Recently, ofF-nuclear tidal flares have also been 
suggeste d in the literature to be probes of recoiling 
SMBHs jKomossa fc MerrittI l2008[ iMerritt et a.1 



20091 : IStone fc Loebl 120111120121 : iLi et al.l l2oT^ 

However, an off-center tidal flare can also be pro- 
duced by SMBHs embedded in merging galaxies. 
In this paper, we calculated the tidal flare rates 
produced by dual SMBHs in a particular evolu- 
tionary stage when the two SMBHs are still un- 
bounded to each other and isolated in the cores 
of merging galaxies. We considered three major 
processes responsible for the loss-cone feeding in 



the merger system, namely, two-body stellar relax- 
ation, tidal perturbation by the companion galaxy, 
and chaotic stellar orbits in triaxial gravitational 
potential. 

By employing an analytical model to calculate 
the stellar disruption rates for both SMBHs in 
the two merging galaxies, we found that prior 
to the formation of SMBHB the stellar disrup- 
tion rate would be enhanced by as large as two 
orders of magnitude in both galaxies. The en- 
hancement is dominated by tidal perturbation and 
occurs when the two galaxies are so close that 
the stars inside the influence radius of the cen- 
tral SMBH are significantly perturbed. We have 
shown that the enhanced stellar disruption rate 
depends on the masses, mass ratio, and density 
profiles of the two galaxies, as well as the distance 
d between the two galaxy cores. In the fiducial 
model with {M7,q,j,l3) = (1, 0.3, 1.75, 2), the en- 
hancement starts when the perturber galaxy ap- 
proaches approximately twice the effective radius 
of the central galaxy {d ~ 2re). In more massive 
systems with M7 > 10, where the stellar disrup- 
tion rates due to two-body relaxation are gener- 
ally lower, the enhancement starts as soon as d 
shrinks to lOrg. As a result, the phase with en- 
hanced stellar-disruption rate extends to an evo- 
lutionary stage much earlier than the formation 
of bound SMBHB, which considerably increases 
the detection rate of wide SMBH pairs in tidal- 
flare surveys. When d shrinks to about the in- 
fluence radius of the central SMBH {d ^ 2ri,), 
the stellar disruption rate in the fiducial model 
increases to 3 x 10^"^ yr^^ in the main galaxy 
and to 2 X 10""^ yr~^ in the satellite. Compared 
to the peak rates in the later evolutionary stages 
with gravitatioiially-bound binary SMBHs (e.g . 
Chen et al.l[2009l 1201 ll: IWegg fc Nate Boddl201lli . 



the total stellar-disruption rates before SMBHs 
become bounded are smaller by only a factor of 
a few. In more massive or equal-mass (q > 0.3) 
mergers, the stellar disruption rates could be even 
higher. 

The above results showed that the tidal disrup- 
tion rates by off-center SMBH pairs in merging 
galaxies are several order magnitude s higher than 
those by recoiling off - nuclear SMBH s ( Komossa fc Merritt 



20081 : iLi et al.|[20ll IStone fc Loebl |2012|) . implv^ 



ing that off-center tidal disruption flares would 
be overwhelmed by the SMBH pairs in merging 
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galaxies. Therefore, it would be challenging to 
distinguish recoiling SMBHs in off-center tidal 
disruption flares. One possible way to distinguish 
the two kinds of off-center tidal disruptions may 
be to identify the evolutionary stages of galaxies. 
Recoiling SMBHs are in galaxies at late stages 
of mergers, while un-bounded SMBH pairs are 
in galaxies at early or middle stages of mergers. 
Early stages of major mergers when galaxies are 
widely separated may be identified with the dis- 
turbed morphology of host galaxies. However, 
morphological signatures of galaxy merger are 
weak during the middle or late stages of major 
mergers, as well as during the whole stages of 
minor mergers, therefore it would be also a chal- 
lenge to identify these merger stages. Another 
difference may be among the properties of star 
clusters around the off-center SMBHs. A recoiling 
SMBH is expected to reside in an ultra-compact 
bound star cluster of mass much smaller than the 
black-hole mass, of size much smaller than the 
black-hole influence radius, and of stellar-velocity 
dispersion much larger than that of host galaxy 
(|Merritt et all 120091 : TTi et al.ll2012[ ). It may also 
associate with a massive cloud of unbound stars, 
whose mass is comparable to the black-hole mass, 
size comparable to the black-hole influence ra- 
dius, and stellar-velocity dispersion comparable 
to or greater th an that of the host galactic nuclei 
(|Li et al.ll2012l ). While the star clusters hosting 
the the secondary black holes in minor mergers 
are the remnants of the tidally truncated satellite 
galaxies. These clusters are orders of magnitude 
heavier than the secondary SMBHs, their sizes are 
much larger than the influence radii of the sec- 
ondary or the primary SMBHs, and their stellar- 
velocity dispersions are comparable to those of 
typical dwarf galaxies but significantly smaller 
than those of the primary galactic nuclei. There- 
fore, the two types of star clusters should differ 
significantly in their sizes, stellar-velocity disper- 
sions, and the mass ratios between SMBHs and 
star clusters, which could be identified with deep 
photomctrical and spectroscopical observations. 

When a pair of SMBHs evolve to about the 
influence radius, d ^ r\y, the enhanced stellar- 
disruption rates can be as high as 10~^ yr~^. For 
such a high tidal disruption rate, multiple tidal 
flares may occur in the same galaxy within a time 
span of decades. Unlike the recurring tidal flares 



in binary or recoiling SMBH systems, the flares in 
merging galaxies are contributed by wide SMBH 
pairs separated by r;, 1 — 10 pc (depending on 
black hole mass and galaxy density profile). Note 
that a separation of 10 pc at redshift z = 0.1 
(1) corresponds to an angular size of 5 (1) mil- 
liarcsec (mas). As a result, spatial offsets be- 
tween successive tidal flares in such a merging sys- 
tem may be detected by instruments such as Gaia 
and Figures [3] and U imply that such flip- 

flop flares could occur in the galaxy mergers with 
10^ Mq < M, < 10^ Mq and q > 0.1. The 
mergers with A/, < 10^ could not produce re- 
curring flares because the stellar-disruption rate 
is too low. When lO"^ Mq < M, < 10^ Mq but 
q < 0.1, most flares are produced in the main 
galaxy; therefore, the recurring flares are unlikely 
to display spatial offset, and would be indistin- 
guishable from those in binary or newly-formed 
recoiling SMBH systems. When M, > 10® Mq, 
the SMBH in the main galaxy would directly cap- 
ture stars, mostly without producing flares, while 
the SMBH in the satellite could still produce tidal 
flares if M,_s < 10® Mq. In the last case, al- 
though the recurring flares occur at the same sky 
position, they should be displaced from the cen- 
ter of the minor merger by an amount of rf, > 10 
pc. Sources with such high flaring rate and large 
off-center displacement cannot be produced by bi- 
nary or recoiling SMBHs. The above discussions 
suggest that with the aid of telescopes with high 
spatial resolution, the cause of the recurring tidal 
flares can be distinguished. 

In the universe, the fraction of tidal flares con- 
tributed by galaxy mergers is proportional to the 
total number of the disrupted stars during merger. 
During each merger, the number of tidal flares 
contributed by phase I, when the separations of 
galaxies are about r^, < d < lOr^, is about 
ni - 103g-i(d/103r6)'3/2 [Figures Q and (0)]. 
The scaling ni (x q~^ implies that ni is determined 
mainly by minor mergers. Suppose a galaxy expe- 
riences N mergers during a Hubble time (^ 10^" 
yr), then during one duty cycle, the number of 
tidal flares contributed by the isolation phase is 
about Us ^ 2 X 10^ /N, if two-body scattering is 
the dominant relaxation process. Since a galaxy 



^See Section 3.7 of LSST Science Book Version 2.0, 
|http: / /www.lsst.org/lsst/scibook| 
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with My = 1 (10) at rcdshift z = has experi- 
enced typicahy N ^ 1 (100) merRers and most 
mergers have g ^ 0.1 ( Hopkins et al. 2010f ). ac- 
cording to the ratio ni : Ug, we find that about 
~ 5% (50%) tidal flares are contributed by phase- 
I galaxy mergers (d ^ lO^rf,). 

For typical mergers with g <C 1, according to 
Figure [SJ the majority of the tidal flares are pro- 
duced in satellite galaxies during phase I, unless 
the satellite galaxies have low surface brightness. 
This result implies that a large fraction of gen- 
uine tidal flares would be displaced by several 
from the centers of the merging systems. Given 
that an offset of 2re ~ 500 pc corresponds to 250 
(60) mas at z = 0.1 (1), these offset tidal flares 
could be misidcntified as supernovae or gamma- 
ray bursts by careless classification schemes. They 
may also be mistaken as "wafceif ' recoiling quasars 
(e.g.. iKomossa fc Merritt 2008 ) or ^^orphan tran- 
sients^'' (X-/7-ray transients either uncorrelated 
with bursts in low-energy bands or without de- 
tectio n of optical counterparts, e.g., iHoran et al 



20091 ) because of the relative dimness of the satel- 
lites. The mis-identification could be very com- 
mon in massive galaxies, because the physical scale 
of re is larger. 

During phase II when the stellar disruption 
rates are enhanced by galaxy mergers, the main 
galaxies would contribution typically more than 
half of the tidal fiares, unless g ~ 1 or the surface 
brightnesses of the main galaxies are low. This re- 
sult indicates that in an advanced merger, where 
the separation between the two galaxy cores is less 
than the effective radius of the main galaxy, the 
tidal flares preferentially reside in the massive nu- 
cleus of the system. According to Figures [7] and 
m the number of tidal fiares contributed by such 
advanced merger phase does not depend on d/rf,, 
and scales as nn ^ 200 (j™Mf, where m < and 
n > are analytical indices derived in § H) There- 
fore, the biggest contribution is expected to come 
from minor mergers in massive systems. Because 
nil is typically less than rij, the contribution of 
tidal fiares from phase II is typically smaller than 
that from phase I. However, for the most massive 
systems with My > 10 in which the main SMBHs 
mostly swallow the stars without producing tidal 
flares, one major merger [q > 0.3), or one mi- 
nor merger {q < 0.3) between galaxies of low sur- 
face brightness, would produce more tidal flares 



in phase II than in phase I. In these particular 
systems, a greater fraction of tidal fiares would be 
contributed by close SMBH pairs with separations 
10 < d/n < 100. 

It is important to know the relative contri- 
butions of tidal flares by single (us), binary 
(rib), and merging SMBH systems (ni and nn). 
The total number of flares produced by a re- 
coiling black holes i s typically smaller than 10' ^ 



( Komossa fc MerrittI 120081: IStone fc Loebl 12011 
Li et al.l bOliT theref ore negligib l e in t he com- 



parison. According to IChen et al.l (|201lh . during 
the lifetime of an SMBHB with g < 1, the in- 
teraction between the binary and the surround- 
ing dense stellar cusp will produce a number 



7 X 



104 ^(2-7)/(6-27)^/i 



11/12 



of tidal 



of rib 

flares. Suppose a galaxy on average experi- 
ences N mergers and M (M < N) of them re- 
sult in the formation of SMBHBs. Then being 
averaged by one duty cycle of galaxy merger, 
TT-s : rif, : (nj -|- nn) is about 20 : 5M : N, 
where we used q = 0.1 because minor merg- 
ers are the m o st common ( Volonteri et al.l 2003 ; 



Hopkins et all l2010t iMcWilhams et al.l boTibh 



tvoical N are fl. 10. 10^) (Hopkins et al. 


2O10l) or 


significantly higher ( McW 


Uiams et al.l 


2013allbl: 


Bedorf & Porteeies Zwart 1 


20131), while M are 



predi cted to be greater than 1 (jVolonteri et al 



20031 ). These numbers highlight the significant 
contribution of tidal flares from merging systems 
with multiple SMBHs. To give more accurate cal- 
culations, one has to combine the cosmic merger 
history of galaxies, as well as the formation rate 
of SMBHBs of different masses and mass ratios. 
Such calculations and the assessment of their un- 
certainties are beyond the scope of the current 
paper and will be addressed in a future paper. 

We arc grateful to Shuo Li, Zuhui Fan, Rainer 
Spurzem, and Thijs Kouwenhoven for helpful com- 
ments. We also thank Alberto Sesana for earlier 
discussions on this topic. This work is supported 
by the National Natural Science Foundation of 
China (NSFC11073002). F.K.L. also thanks the 
support from the Research Fund for the Doctoral 
Program of Higher Education (RFDP), and X.C. 
acknowledges the support from China Postdoc Sci- 
ence Foundation (2011M500001). 



15 



REFERENCES 

Bahcall J. N., & Wolf R. A., 1976, ApJ, 209, 214 

Bedorf, J Portegies Zwar t, S. 2013, MNRAS, in 
press (larXiv:1301.6784p 

Begelman, M. C, Blandford, R. D., & Rees, M. 
J., 1980, Nature, 287, 307 

Binney, J., & Tremaine, S. 2008, Galaetie Dy- 
namics: Seeond Edition, Princeton University 
Press, Princeton, NJ USA 

Bloom, J. S., Giannios, D., Metzger, B. D., et al. 
2011, Science, 333, 203 

Boylan-Kolchin, M., & Ma, C.-P. 2007, MNRAS, 
374, 1227 

Brockamp, M., Baumgardt, H., & Kroupa, P. 

2011, MNRAS, 418, 1308 

Burrows, D. N., Kennea, J. A., Ghisellini, G., et 
al. 2011, Nature, 476, 421 

Cenko, S. B., Krimm, H. A., Horesh, A., et al. 

2012, ApJ, 753, 77 

CentrcUa, J., Baker, J. G., Kelly, B. J., & van 
Meter, J. R. 2010, Rev. Mod. Phys., 82, 3069 

Chen, X., Liu, F. K., & Magorrian, J., 2008, ApJ, 
676, 54 

Chen, X., Madau, P., Sesana, A., & Liu, F. K. 
2009, ApJ, 697, L149 

Chen, X., Sesana, A., Madau, P., & Liu, F. K. 
2011, ApJ, 729, 13 

Cohn, H., & Kulsrud, R. M. 1978, ApJ, 226, 1087 

Colpi, M. & Dotti, M., 2011, Advanced Science 
Letters, 4, 181 

Donley, J. L., Brandt, W. N., Eracleous, M., & 
Boiler, T. 2002, AJ, 124, 1308 

Faber, S. M., Tremaine, S., Ajhar, E. A., et al. 
1997, AJ, 114, 1771 

Frank, J., & Rees, M. J. 1976, MNRAS, 176, 633 

Hills, J. 1975, ApJ, 254, 295 

Hopkins, P. F., et al. 2010, ApJ, 724, 915 



Hopman, C, & Alexander, T., 2006, ApJ, 
645,1152 

Horan, D., Acciari, V. A., Bradbury, S. M., et al. 
2009, ApJ, 695, 596 

Hughes, S. A. 2009, ARA&A, 47, 107 

Gezari, S., Strubbe, L., Bloom, J. S., et al., 2009, 
Astro2010: The Astronomy and Astrophysics 
Decadal Survey, Science White Papers, no. 88 

Ivanov, P. B . & Chernyakova, M. A., 2006, A&A, 
448, 843 

Ivanov, P. B., Polnarev, A. G., & Saha, P., 2005, 
MNRAS, 358, 1361 

Jiang, C. Y., Jing, Y. P., Faltenbacher, A., Lin, 
W. P., & Li, C. 2008, ApJ, 675, 1095 

Kesden, M. 2012, Phys. Rev. D, 85, 024037 

Khan, F., Just, A., & Merritt, D. 2011, ApJ, 732, 
89 

Komossa, S., 2002, Rev. Mod. Astron., 15, 27 

Komossa, S., & Merritt, D., 2008, ApJ, 683, 21 

Lauer, T. R., et al. 2005, AJ, 129, 2138 

Li, S., Liu, F. K., Berczik, P., Chen, X., & 
Spurzem, R. 2012, ApJ, 748, 65 

Lightman, A. P. & Shapiro, S. L., 1977, ApJ, 211, 
244 

Liu, F. K., Li, S., & Chen, X., 2009, ApJ, 706, 
L133 

Magorrian, J. & Tremaine, S., 1999, MNRAS, 309, 
447 

Marconi, A., & Hunt, L., 2003, ApJ, 589, 21 

McWiUiams, S.T., Ostriker, J.P., & Pretorius, F. 
2013a, Science, submitted (larXiv:1211.4590p 

McWiUiams, S.T., Ostriker, J. P., & Pretorius, F. 
2013b, ApJ, submitted (|arXiv:1211.5377p 

Merritt, D. 2006, ApJ, 648, 976 

Merritt, D., & Cruz, F. 2001, ApJ, 551, L41 

Merritt, D. & Milosavljevic, M., 2005, Liv. Rev. 
Rel., 8, 8 



16 



Merritt, D. & Poon, M. Y., 2004, ApJ, 606, 788 

Merritt, D., Schnittman, J. D., & Komossa, S. 
2009, ApJ, 699,1690 

Merritt, D., & Vasiliev, E. 2011, ApJ, 726, 61 

Ostrikcr, J. P., Spitzer, L., & Chavalier, R. A., 
1972, ApJ, 176, L47 

Perets, H. B., Hopman, C, & Alexander, T., 2007, 
ApJ, 656, 709 

Poon, M. Y., & Merritt, D., 2001, ApJ, 549, 192 

Poon, M. Y., & Merritt, D., 2004, ApJ, 606, 774 

Preto, M., Berentzcn, I., Berczik, P., & Spurzcm, 
R. 2011, ApJ, 732, L26 

Ranch, K. P., & Ingalls, B., 1998, MNRAS, 299, 
1231 

Ranch, K. P., & Trcmaine, S. 1996, NewA, 1, 149 

Rees, M. J. 1988, Natnre, 333, 523 

Roos, N., 1981, A&A, 104, 218 

Schnittman, J. D. 2010, ApJ, 724, 39 

Seto, N., & Muto, T. 2011, MNRAS, 415, 3824 

Stone, N., & Loeb, A. 2011, MNRAS, 412, 75 

Stone, N., & Loeb, A. 2012, MNRAS, 422, 1933 

Syer, D. & Ulmer, A., 1999, MNRAS, 306, 35 

Tremainc, S., et al. 2002, ApJ, 574, 740 

van Velzen, S., Farrar, G. R., Gezari, S., et al. 
2011, ApJ, 741, 73 

Volonteri, M., Haardt, F., & Madau, P. 2003, ApJ, 
582, 559 

Wang, J.-X. & Merritt, D., 2004, ApJ, 600, 149 
Wegg, C, & Nate Bode, J. 2011, ApJ, 738, L8 
Young, P. J. 1977, ApJ, 215, 36 



This 2-column preprint was prepared with the AAS lAT^X 
maeros v5.2. 

17 



